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ABSTRACT 

We have carried out photometric follow-up observations of bright transiting extrasolar planets using 
the CbNUOJ 0.6m telescope. We have tested the possibility of obtaining high photometric precision 
by applying the telescope defocus technique allowing the use of several hundred seconds in exposure 
time for a single measurement. We demonstrate that this technique is capable of obtaining a root- 
mean-square scatter of order sub-millimagnitude over several hours for a V ~ 10 host star typical for 
transiting planets detected from ground-based survey facilities. We compare our results with transit 
observations with the telescope operated in in-focus mode. High photometric precision is obtained 
due to the collection of a larger amount of photons resulting in a higher signal compared to other 
random and systematic noise sources. Accurate telescope tracking is likely to further contribute to 
lowering systematic noise by probing the same pixels on the CCD. Furthermore, a longer exposure time 
helps reducing the effect of scintillation noise which otherwise has a significant effect for small-aperture 
telescopes operated in in-focus mode. Finally we present the results of modelling four light-curves for 
which a root-mean-square scatter of 0.70 to 2.3 milli-magnitudes have been achieved. 

Key words : extrasolar planets, transiting planets, data modelling, photometric noise, defocus 
technique 


1. Introduction 

The first discovery of an extrasolar planet in 1995 
(Mayor & Queloz 1995) resulted in the opening of a 
completely new astronomical research area. However, 
detailed information of the planet obtained via radial 
velocity measurements is limited as only the minimum 
mass (among other parameters) could be inferred due 
to an unknown orbital inclination. The class of tran¬ 
siting extrasolar planets has changed this ambiguity. 
The discovery of a transiting extrasolar planet (TEP) 
(Charbonneau et al. 2000; Henry et al. 2000) allows 
the determination of physical properties of the planet 
(mass, radius among others) and its host star (limb- 
darkening, effective temperature among others) when 
combined with spectroscopic observations. In princi¬ 
ple, obtaining precise measurements of these properties 
will distinguish a given TEP to be a rocky terrestrial- 
size or gaseous Jupiter-size planet by inferring its ab¬ 
solute size and mass. Therefore accurate properties 
will help to constrain planet formation theories. One 
possible technique to obtain high-precision photometric 
light-curve is the use of telescope defocus technique and 
has been successfully applied in various studies South- 
worth et al. (2009). This research presents results from 
our attempt to obtain precise light-curves of TEPs us- 
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ing a 0.6m telescope. In particular we present results 
for which the telescope was operated in-focus as well as 
out-of-focus. In section 2 we outline the background of 
the defocus technique. In section 3 we present our ob¬ 
servations of seven TEPs observed since early 2012. A 
description of data reduction and light-curve modelling 
is given in section 4 and 5. Our results and analysis are 
given in section 6 and 7 followed by a conclusion in sec¬ 
tion 8. 

2. Defocus photometry 

Traditional astronomical photometry aims to mea¬ 
sure the brightness of a star with the telescope kept in 
focus resulting in a well-defined point-spread function 
(PSF). This is usually desired when structural infor¬ 
mation is needed (e.g. star cluster observations to re¬ 
solve individual stars) or when the observed star field 
is crowded to avoid confusion due to overlapping PSFs. 
However, in certain situations the operation of the tele¬ 
scope in a defocus mode opens up for the possibility to 
achieve an increased signal-to-noise (S/N) ratio result¬ 
ing in a higher photometric precision. Nevertheless the 
same noise processes are in operation whether the tele¬ 
scope is in-focus or de-focused. To obtain high-S/N 
measurements over an extended time period both in¬ 
strumental and environmental effects are important to 
consider. 
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In general several types of noise sources contributes 
to the overall error budget of a single measurement. 
The first class is random noise (also known as statis¬ 
tical, white or time-uncorrelated noise) and often re¬ 
ferred to Gaussian noise associated with Poisson count¬ 
ing statistics of photons. CCD read-out noise, CCD 
dark-current, shot-noise of the background as well as 
the target star are all referred to random noise. In 
a special case atmospheric scintillation noise is also of 
random nature, but depending on other factors such as 
telescope aperture (see below for more details). To de¬ 
crease white noise it requires an increase in the count¬ 
ing statistic of photons. The second class of noise 
source is of systematic nature and harder to identify 
and/or quantify. Examples of systematics would be 
selecting aperture sizes for photometry or flat-fielding 
errors involving the positioning of the target star on 
the same pixels. Controlling the latter is ultimately 
linked with telescope pointing stability and involves 
how well the telescope is tracking over an extended 
time period. Also the choice of comparison stars for 
differential photometry is a potential source of system¬ 
atic noise. Furthermore, the telescope optical align¬ 
ment (astigmatism) and general telescope optics also 
are potential sources for introducing systematic noise. 

Other noise sources include non-random noise (also 
known as red, 1// or time-correlated noise) often (but 
not exclusively) associated with atmospheric extinction 
encountered at the beginning and/or end of a transit 
observation when the target star is observed through 
high air-nrass. The final error source is related to the 
intrinsic variability of the star itself such as star spots, 
flares and pulsations. All of these noise sources have 
the effect of deteriorating the photometric time-series 
measurements in one or the other way making it dif¬ 
ficult to obtain a photometric precision of less than a 
rnilli-magnitude. 

One key aspect of limiting the S/N of a photometric 
measurement is the dynamic range of a CCD detec¬ 
tor dictating the peak flux and thereby setting a nat¬ 
ural upper limit on the exposure time for a telescope 
operated in a in-focus mode. CCD detectors with a 
higher dynamic range allow for a longer exposure time 
thereby circumventing the danger of reaching the sat¬ 
uration limit. Ideally, in order to obtain a high signal- 
to-noise (S/N) ratio plenty of photons are needed (oth¬ 
erwise achieved by using a large mirror aperture). The 
larger the photon counts the lower the Poisson noise. 
One practical way to allow for a longer exposure time 
without saturating the CCD is to spread the light over 
many pixels by defocusing the telescope. As a result 
the stellar PSF of the target (and all other stars in the 
same field) broadens and decreases the peak value of 
an exposed pixel. 

In Fig. 1 we show an example of two stars of 
fairly similar brightness, but defocused for the case of 
HATP22. Simply, in a defocus mode, the number of 
photons are distributed over a larger area of the CCD 
allowing a prolonged exposure and hence a stronger 




Fig. 1 .— Surface plots of two PSFs. Upper panel: PSF 
of X03 (2012/02/11) with the telescope well focused and 
an exposure time of 12s. Lower panel: PSF of HATP22 
(2014/01/13) with the telescope defocused heavily allowing 
for a much longer exposure time of 195s. Both stars are of 
same brightness (see Table 2). 


signal to build up and stand out from the baseline of 
combined noise sources. 

However, a natural upper limit of exposure time 
{t eX p ) is determined by the single frame read-out time 
( t ro ) and the transit duration ( td U r ) as follows 

Tdur — N * (texp T t ro ) (1) 

where N is the number of frames recorded during the 
transit. To characterise a planetary transit from a de¬ 
tailed modeling process the light-curve sampling has to 
be sufficiently high in order to ensure enough measure¬ 
ments. As a rule of thumb we aim to obtain ./V = 30 to 
40 measurements during a given transit event resulting 
in 4 to 6 measurements of the ingress and egress phase. 

Several advantages and disadvantages of the defocus 
technique contribute to either decrease or increase the 
S/N ratio. The variable atmospheric seeing effect and 
local temperature changes are less important for obser¬ 
vations with a broad PSF. However, a high atmospheric 
transparency is still desireable. Distributing the pho¬ 
tons over many pixels also has the effect of decreasing 
the photometric noise due to intrinsic CCD pixel-to- 
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Fig. 2. — Expected photometric error due to scintillation 
noise as a function of exposure time and mirror aperture 
valid for a telescope at 80m above sealevel. Top panel: Con¬ 
sidering an observation through air-mass 1.0. Lower panel: 
Considering air-mass of 1.5. See text for more details. 


pixel response variations. The total Hatfield noise con¬ 
tribution within a larger signal aperture (contribution 
of each pixel added in quadrature) is lower when com¬ 
pared to a similar in-focus measurement. In essence, 
the flatheld error of the final photometric measurement 
should average out using a larger number of indepen¬ 
dent pixels due to its random nature. One technical 
possibility along which one can increase the number of 
pixel is to operate the CCD in one-by-one binning. As 
the maximum number of pixels are requested the read¬ 
out time in single-pixel binning mode is maximised. 
Forming super pixels via two-by-two (or higher) bin¬ 
ning would decrease the read-out time. However, a dis¬ 
advantage of defocus photometry is the contribution of 
CCD read-out noise. The more pixels are exposed the 
larger the total read-out noise on a given measurement. 
In addition, the contribution of background noise in the 
broadened PSF also will decrease the photometric S/N. 

A positive effect of a prolonged exposure for a small- 


aperture telescope is the beating of scintillation noise. 
Scintillation is caused due to the refraction of starlight 
(twinkle) by turbulent cells in the atmosphere. The ef¬ 
fect introduces a time-varying fluctuation in the star’s 
brightness. Its effect is largest for telescopes with small 
mirror diameters and short exposure times. Therefore, 
scintillation will set a lower-limit of the photometric 
precision obtained for a small-aperture telescope. The 
first to model the effect of scintillation for bright stars 
is Young (1967). Using the expression given in Hart¬ 
man et al. (2005) we have examined the dependence of 
scintillation noise on telescope aperture, air-mass, tele¬ 
scope height and exposure time. In Fig. 2 we plot the 
scintillation noise (in milli-magnitude) as a function of 
exposure time for three different telescope apertures 
located at 80m above sealevel. The top panel shows 
the relationship for air-nrass=1.0 and the lower panel 
for air-mass=1.5. In general, the scintillation noise de¬ 
creases for either increasing telescope aperture or for 
increasing exposure time at a given air-nrass. Consid¬ 
ering a 60cm telescope, at unit air-mass and integrat¬ 
ing for 250 seconds the scatter per datapoint expected 
from scintillation is below 0.3 mnrag. At air-mass 1.5 
this limit increases to 0.6 nnnag - more than sufficient 
to record transits depths of 1 - 2% with good precision. 
However, scintillation imposes a lower limit which is 
achievable for only bright stars. 

3. Observations 

All observations* presented in this work were carried 
out using the 0.6m telescope (hereafter CbNUOJ) lo¬ 
cated in Jincheon at an altitude of 87m above sealevel. 
The telescope was installed by Korea Astronomy and 
Space Science Institute (KASI) and is operated by 
Chungbuk National University Observatory, Republic 
of Korea. The telescope optics follows the Richey- 
Chretien design attached to an equatorial mount. Tele¬ 
scope tracking is performed via a servo motor control. 
More details of this telescope was described by Kim 
et al. (2014). 

During the observing period the telescope was equip¬ 
ped with two different CCD imagers located at the 
Cassegrain f/2.92 prime focus and installed by the Ko¬ 
rea Astronomy and Space Science Institute (KASI). 
Initially we used the 1530 x 1020 pixel SBIG ST-8XE^ 
which was later replaced (currently installed) by the 
more advanced 4096 x 4096 pixel SBIG STX-16083* 
CCD camera. The former camera had a field of view 
of 27 x 18 and the current camera provides a field of 
view of 72 x 72 . Both cameras had a pixel scale of 1.05 
arcsec/pixel. Single Johnson/Cousins UBVRI filters^ 


* photometric data can be obtained in ASCII form from the cor¬ 
responding author 

t http: / / archive.sbig.com / sbwhtmls / st8.htm 
^ https: / / www.sbig.com / products / cameras / stx/stx-16803 
§ www.astrodon.com 
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Table 1. 

Observation log of seven TEPs. The transit of HATP09 was the first transit to be recorded using the new ST-16083 

camera. 


Target 

Date 

Start time 

End time 

Exposure time 

Filter 

# Obs. 

Air mass 

Camera 

XQ5 

2012 /02/10 

UT 10:09 

UT 18:25 

75s - 105s 

R 

288 

1.00-1.56 

ST-8 

X03 

2012 /02/11 

UT 09:58 

UT 15:06 

10 s - 12s 

R 

1084 

1.08-1.86 

ST-8 

X04 

2012/02/17 

UT 10:55 

UT 18:05 

27s - 35s 

R 

568 

1.12-1.63 

ST-8 

HATP25 

2012/11/15 

UT 10:49 

UT 18:50 

170s - 190s 

R 

149 

1.02-1.83 

ST-8 

HATP09 

2012/12/23 

UT 12:33 

UT 20:22 

230s (fixed) 

R 

109 

1.00-1.62 

ST-16083 

HATP22 

2014/01/13 

UT 14:41 

UT 21:23 

195s (fixed) 

R 

114 

1.03-1.35 

ST-16083 

HATP12 

2014/03/09 

UT 13:42 

UT 18:39 

260s (fixed) 

R 

65 

1 .00-1.66 

ST-16083 


Table 2. 

Details of seven host stars and their known TEPs. Data was obtained from http://simbad.u-strasbg.fr/simbad/sim-fid and 

http: / / exoplanet, eu. 


Target 

RA (J2000) 

Dec (J2000) 

V (mag) 

Depth (mag) 

Duration (min) 

X05 

07 46 51.95 

+39 05 40.5 

12.1 

0.014 

193 

X03 

04 21 52.71 

+57 49 01.9 

9.9 

0.0048 

173 

X04 

07 21 33.17 

+58 16 05.0 

10.8 

0.0108 

264 

HATP25 

07 21 33.17 

+58 16 05.0 

13.2 

0.0204 

169 

HATP09 

07 20 40.48 

+37 08 26.5 

12.3 

0.0126 

206 

HATP22 

10 22 43.59 

+50 07 42.0 

9.8 

0.0119 

172 

HATP12 

13 57 33.48 

+43 29 36.7 

12.8 

0.0204 

140 


are available via an electronic filter wheel system. For 
most of our observations we used the Cousins R-filter. 
In some observations we experimented with clear fil¬ 
ters. The large field of view in the current instrumental 
setup is ideal for differential photometry when carrying 
out ensemble photometry analysis (Everett & Howell 
2001). Both cameras are cooled electronically. 

Unbinned CCD read-out time for the ST-8XE was 
about 10 seconds while for the STX-16083 camera it 
was around 16 seconds. The two relatively short read¬ 
out times are favourable for high-cadence transit pho¬ 
tometry and allowed us to stay on target for a larger 
fraction of time for photon collection. Another benefit 
of short read-out times is a smaller contribution to the 
noise budget due to read-out noise imposed on a sin¬ 
gle frame as the high photon counts will dominate for 
defocus measurements. In the beginning (2012 season 
using the ST-8 camera) we followed a strategy where 
we changed the exposure time during a given night ac¬ 
cording to weather conditions. Later we abandoned 
this strategy and kept the exposure time fixed for con¬ 
sistency of the recorded data. Usually to our experi¬ 
ence even thin clouds will have a temporal dramatic 
effect on the sampled light-curve when trying to aim 
for a mean ~ 1 mrnag photometric precision over 3 
to 4 hours. Changing exposure times by a few percent 
will not have a significant impact on the signal-to-noise 
ratio and our goal was to ensure a well-sampled light- 


curve. All observations were carried out with the cam¬ 
eras operated in unbinned mode providing a maximum 
number of pixels favourable for defocus photometry. 
Prior to the defocus observations we aimed to obtain 
2-3 in-focus images in order to identify any nearby com¬ 
panions. The average seeing at CbNUOJ is around 2 
arcsec. 

Calibration frames (bias and skyflats) were obtained 
for each night using appropriate filter. We carried out 
tests for which the science images were either calibrated 
or kept in their original form. We will report about 
the results from those experiment in the next section. 
An observation log is given in Table 1. Details of the 
observed host stars are given in Table 2. 

4. Data reduction 

We used the DEFOT^ data reduction pipeline outlined 
in Soutliworth et al. (2009). The software is written in 
IDL and makes use of the NASA IDL astronomy library 
package astrolib. The method of obtaining photo¬ 
metric measurements is following the algorithm de¬ 
scribed in the DA0PH0T photometry package. Aperture 
photometry is done using the astrolib/aper routine. 
UTC (Coordinated Universal Time) time stamps in 
the format YYYY-MM-DD/HH-MM-SS are obtained 


■ htt p: //www. astro. keele. ac. uk/~j kt 
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from each image FITS header. Actual times for a 
single photometric measurement were taken to be the 
mid-exposure time. To convert to Barycentric Julian 
Date (BJD) in the Barycentric Dynamical Time (TDB) 
standard we used the utc2bj d. pro conversion routine^ 
Eastman et al. (2010). BJD time stamps expressed in 
the TDB standard will be particularly useful in when 
searching for transit timing variations (TTV). 

Apertures for the inner and background annulus 
were fixed (but their optimum values determined from 
experimentation) in all frames for each data set. We 
noticed a significant frame-to-franre shift due to track¬ 
ing imperfections of the telescope. To ensure that the 
apertures are following the selected stars (target and 
comparison) we made use of the image registration 
cross-correlation algorithm available in the DEFOT pack¬ 
age. For each program image a summed pixel row and 
pixel column intensity profile is calculated and com¬ 
pared with the same profile of a reference frame. A 
least-squares high-order polynomium fit is then per¬ 
formed to calculate any shifts and a given frame is then 
corrected accordingly in order to match the reference 
frame. The downside of the package is that only a sin¬ 
gle aperture for the three annuli is applied to all stars. 
In practice this limits the selection of comparison stars 
similar in brightness (and hence angular extend on the 
CCD) with the target. From experimentation we used 
apertures that yielded the smallest root-mean-square 
(RMS) scatter around a best-fit model (see next sec¬ 
tion). Table 3 gives an overview of the various aper¬ 
tures used for each target. 

Comparison stars were checked for internal variabil¬ 
ity and none were found for the selected stars. The 
DEFOT package allows to carry out ensemble differential 
photometry using weighted flux summation. This pro¬ 
cedure minimises the overall Poisson noise while also 
avoids distortion in the transit light curve. More de¬ 
tails of this procedure can be found in Southworth et al. 
(2009) and Everett & Howell (2001). 

Long-term trends over the observing window were 
removed using a polynomium fit to out-of-transit data. 
The resulting light curve was then normalized to zero 
differential magnitude. In all normalisations we applied 
a first-order polynomium fit. Higher-order terms would 
have the potential of introducing artifical distortions to 
the transit light curve. 

In some cases we noticed obvious data outliers most 
likely due to passing thin clouds. These were removed 
prior to our photometric analysis. To find stable com¬ 
parison stars we selected all suitable stars and carried 
out a first selection by eye. Then a subset of stable stars 
(4-15) were included in the final photometric analysis. 
To find optimum apertures we varied the three radii 
in a systematic way. For each experiment we found 
a best-fit solution with associated RMS scatter of the 
residuals. Radii resulting in a least scatter were cho- 


II http://astroutils.astronomy.ohio-state.edu/time/ 


Table 3. 

Final aperture radii used for photometric measurements. 
Ri is the inner aperture. R 2 , R 3 are the apertures for the 
background measurement. All numbers are in pixels. The 
apertures were determined in order to minimise the RMS 
scatter of the data. 


Target 

Ri 

R 2 

i?3 

HATP25 

9 

13 

20 

HATP09 

12 

40 

45 

HATP22 

25 

35 

45 

HATP12 

10 

30 

45 


sen. These experiments were done using raw science 
frames. Using the most optimal apertures we also cal¬ 
ibrated (bias subtraction and flatfielding) the science 
frames to determine any difference in RMS. In some 
cases we found that a calibration resulted in a slight 
improvement in the photometry in other cases we found 
no improvement. 

5. Light-curve modelling 

Light-curves of all transiting planets were modelled 
using the JKTEB0P_V25** least-squares minimisation 
code Southworth (2008); Southworth et al. (2009). 
The code models the light curve of detached eclips¬ 
ing binary stars using biaxial spheroids. Fundamental 
parameters (star and planet are referred to with sub¬ 
scripts A and b, respectively) describing the light-curve 
of a transiting planets are the orbital period P, the time 
of minimum light To (or ephemeris epoch when mod¬ 
elling multiple light-curves simulteanously), the frac¬ 
tional radii of the star (ta = Ra/o, where a is the 
orbital semi-major axis) and planet (n, = Rb/a), the 
orbital inclination (if,) and limb-darkening coeffiecients. 

In practice, the parameters ta and are parameterised 
as their sum [ta + rb) and ratio (k = Vb/rA)- 

In all modelling work we kept the radii sum, radii 
ratio, orbital inclination, the time of minimum light 
and the scale factor (normalised flux level) freely ad¬ 
justable. Since this work presents transit data of a 
single event for a given planet we did not fit for the or¬ 
bital period and kept it fixed. Stellar limb-darkening is 
implemented in JKTEBQP in form of several parametric 
laws. A linear or quadratic limb-darkening law intro¬ 
duces another one or two parameters. Initial guesses of 
known quantities have been obtained from the respec¬ 
tive discovery papers for each system. 

When using JKTEB0P_V25 in its simplest form (TASK3) 
the parameter uncertainties are obtained from the best- 
fit covariance matrix. Therefore they are formal (and 
probably too optimistic) errors. More proper parame¬ 
ter uncertainties are obtained through the implemen- 


**http://www.astro.keele.ac.uk/~jkt 
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BJD - 245000.0 

Fig. 3. — Light-curve of X05 as observed on 2012/02/10 
with the telescope in focus. Upper panel: The RMS scat¬ 
ter around the best-fit model is 2.98 mmag. The large 
scatter does not allow for a detailed limb-darkening treat¬ 
ment. Ingress and egress phases are not characterised well. 
The mid-transit time is at BJD 2,455, 968.13350 ± 0.00061. 
Lower panel: Various choices of exposure times during the 
observing period. Some correlations between short expo¬ 
sures and large scatter is visible. 


tation of bootstrapping, Monte-Carlo and residual- 
permutation (RP) algorithms. To obtain robust pa¬ 
rameter uncertainties for light-curves with correlated 
red noise the RP algorithm is most suitable and avail¬ 
able via TASK 9. 

6. Results - 1 

Our first few test observations with the telescope in 
focus aimed at recording transit light-curves of X05 
(Burke et al. 2008), X03 (Johns-Krull et al. 2008) and 
X04 (McCullough et al. 2008). The resulting light- 
curves are shown in Fig. 3, Fig. 4 and Fig. 5, respec¬ 
tively. In all three cases the transit is clearly visible. 
The photometric point-to-point scatter is no less than 
~ 3 mmag with maximum exposure times ranging from 
105 seconds for X05 (V = 12.1 mag) to 12 seconds for 
X03 (V = 9.9 mag). Short exposure times result in 
a high sampling cadence. The residuals in each plot 




Fig. 4. — Light-curve of X03 as observed on 2012/02/11 
with the telescope in focus. Upper panel: The RMS scatter 
around the best-fit model is 3.18 mmag. The mid-transit 
time was determined as BJD 2,455, 969.03386 ± 0.00024. 
Lower panel: Air-mass versus time. Due to high air-mass at 
the end of observing night the photometric scatter increases. 


all show some degree of systematic effects indicating 
that correlated (red) noise is the most dominant er¬ 
ror source. In all cases the photometric precision is 
too low for a sufficient characterisation of the tran¬ 
sit. The derived planet parameters have a large un¬ 
certainty. Nevertheless we modelled the lightcurves in 
order to obtain an estimate of the mid-transit time. 
Limb-darkening treatment is not meaningful since no 
such effect is readily apparent in the time-series pho¬ 
tometry. Errorbars per observation has been omitted 
in the figures to highlight some details concerning our 
observing strategy, sky condition and telescope track¬ 
ing properties. 

In Fig. 3 we also plot the used exposure time during 
the observing window. Depending on sky condition we 
varied the exposure between 75 seconds and 110 sec¬ 
onds. After BJD 2,455,968.15 the scatter systemati¬ 
cally shifts upwards in the differential magnitude plot. 
This is coinciding with a period of time where the ex¬ 
posure time was chosen to be minimal demonstrating 
that decrease in exposure time deteriorates the photo¬ 
metric precision. This finding motivated us to stick to 
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Fig. 5. — Light-curve of X04 as observed on 2012/02/17 
with the telescope in focus. Upper panel: The RMS scatter 
around the best-fit model is 3.34 mmag. The mid-transit 
time was determined as BJD 2,455, 975.088±0.002. Middel 
and Lower panels: Telescopic X and Y drift of program 
image relative to reference image. 


a constant exposure time during the course of observa¬ 
tion resulting in a homogeneous data set. 

In Fig. 4 we chose to plot the air-mass during the 
X03 transit. We demonstrate that observing through 
a relatively high air-mass (>1.45) results in a signif¬ 
icant increase in photometric scatter most likely due 
to a larger scintillation noise. To obtain high-precision 
photometry of a transiting planet one should aim to ob¬ 
serve through as low air-mass as possible. This require- 



Fig. 6. — Light-curve of HATP25 as observed on 
2012/11/15 with the telescope de-focused to allow for longer 
exposures. The RMS scatter around the best-fit model is 
2.30 mmag with Xr = 1-01. No calibration frames were 
applied. 



BJD - 2450000.0 


Fig. 7. — Light-curve of HATP09 as observed on 
2012/12/23 with the telescope de-focused to allow for longer 
exposures. The RMS scatter around the best-fit model is 
1.68 mmag with \r = 1-42. A linear limb-darkening law 
was determined to model the data best. No calibration 
were performed. 


rnent would exclude targets located on more southern 
latitudes. 

A last characteristic is concerned with the pointing 
ability of the telescope itself. The cross-correlation al¬ 
gorithm as implemented in DEFQT for image registration 
records the number of pixel with which each program 
image was shifted in order to match the reference im¬ 
age. In Fig. 5 we plot the telescopic x- and y-drift 
measured in pixels during the X04 transit observation. 
The x-drift corresponds to ± righ-ascension while y- 
drift is ± declination. In both directions the maximum 
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Fig. 8. — Light-curve of HATP22 as observed on 

2014/01/13 with the telescope heavily de-focused to allow 
for longer exposures. The RMS scatter around the best-fit 
model is 0.70 mmag with \r = 1.21. Bias subtractions and 
flatfielding calibration were performed. 



Fig. 9. — Light-curve of HATP12 as observed on 

2014/03/09 with the telescope slightly de-focused to allow 
for longer exposures. The RMS scatter around the best-fit 
model is 1.85 mmag with Xr = 0.91. No bias subtractions 
and flatfielding calibration were performed. 


drift is around 10 pixels (10.5 arcseconds). While the 
y-drift appears to be of continous nature the x-drift 
occurs in jumps. Pointing imperfections can either be 
due to slag in the motor drives/friction-disk, inaccu¬ 
rate encoder and/or software control. The ability to 
keep the target on the same pixels during the tran¬ 
sit observations contributes to minimising flatfielding 
errors. As an example the recently commissioned K2 
mission (Kepler space telescope extended mission) will 
have its photometric precision decreased significantly 
(compared to the original Kepler mission) due to de¬ 
creased pointing accuracy. 


Table 4. 

Theoretical limb-darkening coeffiencients for a quadratic 
LD law valid for a Cousin R-filter. See text for details. 


Host star 

linear ( u a ) 

quadratic ( v a ) 

HATP25 

0.3879 

0.2906 

HATP09 

0.5588 


HATP22 

0.4285 

0.2646 

HATP12 

0.5669 

0.1687 


7. Results - 2 

In order to obtain high-precision photometric mea¬ 
surements of a transiting planet we then carried out 
tests with long exposure times per observation frame. 
A significant amount of telescope defocus would be re¬ 
quired for the more brighter targets to ensure ADU 
counts to be in the linear regime of the CCD. We ob¬ 
tained four light-curves of four different systems on 
four different nights. Three of the recorded light- 
curves makes use of the new ST-16083 CCD cam¬ 
era. In Fig. 6 to Fig. 9 we show the light-curves of 
HATP25 (Quinn et al. 2012), HATP9 (Shporer et al. 
2009), HATP22 (Bakos et al. 2011), HATP12 (Hart¬ 
man et al. 2009). Light-curves of HATP12 were previ¬ 
ously recorded several times using the Korean lrn opti¬ 
cal telescope at Mt. Lemmon Optical Astronomy Ob¬ 
servatory (LOAO), Arizona, USA Lee et al. (2012). All 
light-curves were observed through a Cousins R-filter. 
The sampling frequency is now due to a longer exposure 
time. Exposure times of over 3 minutes have been used 
on average. The effect is immanent and we see a signifi¬ 
cant increase in photometric precision compared to our 
first test observations (see previous section). Using a 
fixed exposure time of 260s we achieved a RMS scatter 
of ~ 1.85 mmag for the faintest host star (HATP12, V 
= 12.8 mag). A remarkable RMS scatter of 0.70 mmag 
was obtained for the brightest target HATP22 (V = 
9.8 mag). 

7.1 Treatment of limb-darkening 

Quadratic limb-darkening coefficients for Cousin R- 
filter were obtained using the JKTLDU code which out¬ 
puts theoretically calculated limb-darkening strengths 
for a variety of parametric laws. The coefficients are 
calculated from bilinear interpolation (in effective tem¬ 
perature and surface gravity) from published tables cal¬ 
culated from stellar model atmospheres. In this work 
we obtained the coefficients for the quadratic limb- 
darkening law from tables published in Claret et al. 
(2000). In all cases we assumed the metallicity of 
the host-star to be zero and the atmospheric micro¬ 
variability parameter were set to unity. Using nominal 
published values for T e //,logg, [Fe/H] and v m i C does 


''http: //www.astro.keele.ac.uk/jkt/codes/jktld.html 
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Frequency (cycles/day) 

Fig. 10. — Top panel: Light-curve of HATP16 as ob¬ 
served on the night of 2011/11/03 with the lm reflector 
at Mt. Lemmon Optical Astronomy Observatory (LOAO). 
The RMS scatter around the best-fit model is 2.41 mmag. 
Lower panel: Power spectrum (normalised to the highest 
peak) of the residuals between observations and best-fit 
model. Clearly a 1/frequency dependency is evident. 


not allow one to determine suitable coefficients. 

7.2 Treatment of parameter uncertainties 

JKTEBOP allows parameter uncertainties to be ob¬ 
tained from two different methods depending on the 
noise character in the data set. If the observational 
noise is time-correlated (red noise or covariant noise, 
see eventually Gillon et al. (2009)) then TASK9 is rec¬ 
ommended and makes use of the residual-permutation 
(a variant of prayer-bead) algorithm. If the noise is 
Gaussian or white noise then a standard Monte-Carlo 
algorithm is implemented via TASK8. The former 
method requires the computation of n shifted fits where 
n is the number of data points. The latter method 
requires a larger number of simulations and we used 
10,000 experiments to obtained parameter uncertain¬ 
ties from TASK8. The remaining question is how to 


judge whether a given data set contains red or white 
noise. One method is to calculate the power spectral 
density of the photometric time-series or its residu¬ 
als from a best-fit solution. If the observations have 
a high content of time-correlated (red) noise then the 
power spectrum of the timer-series will exhibit a 1 /fre¬ 
quency dependence: most power at low frequencies. 
To test this behaviour we calculated the power spec¬ 
trum of a single HATP16 transit observed with the 
lm reflector at Mt. Lemmon Optical Astronomy Ob¬ 
servatory. Figure 10 shows the observed transit along 
with a best-fit solution. Systematic trends in the resid¬ 
uals are clearly visible during the observing window. 
The lower panel shows the power spectrum computed 
from a Lomb-Scargle algorithm as implemented in the 
PERI0D04M software package (Lenz & Breger 2005) 
with a clear 1/frequency trend pointing towards a sig¬ 
nificant amount of correlated noise (also visible in the 
residual plot in Fig. 10). We have calculated power 
spectra for the four planetary transits obtained with 
the CbNUOJ telescope and display the results in Fig¬ 
ure 11. In all cases (with the HATP25 transit as a pos¬ 
sible exception with a marginal content of red noise) 
the power spectrum is flat (certainly no 1/frequency 
dependency) and hence no time-correlated red noise is 
present in these data sets. We have therefore decided to 
use JKTEBOPs TASK8 for all transit observations to 
estimate parameter uncertainties. Since HATP25 could 
have red noise components we have also estimated pa¬ 
rameter uncertainties using TASK9. 

For HATP25 we first fitted for the two limb-darkening 
coefficients separately. We find that the data does not 
allow a simulteaneous fitting for the coefficients and 
found that the sum of final coefficients is larger than 
unity (implying an unphysical brightening at the host- 
star limb). We therefore decided to fix those param¬ 
eters at their theoretical values as shown in Table 4. 
The final best-fit parameters for the HATP25 transit 
is shown in Table 5. Since we estimated parameter 
uncertainties using TASK8 and TASK9 we can now 
judge that some amount of time-correlated red noise 
is present for that dataset. The errors obtained from 
TASK9 are consistently larger than for TASK8. 

For HATP09 we also fitted for the two limb-darkening 
coefficients separately, but encountered the same prob¬ 
lem as for HATP25 (i.e unphysical limb-darkening co¬ 
efficients). We then set the two linear and quadratic 
terms to their tabulated values as obtained from JKTLD 
and fixed those parameters. In this case we encoun¬ 
tered that the inclination was kept artificially at 89.96 
degrees to ensure numerical stability of the minimisa¬ 
tion algorithm. However, the best-fit inclination in the 
discovery paper Shporer et al. (2009) found an inclina¬ 
tion of 86.5 degrees. We therefore decided to model 
the light-curve with a linear limb-darkening law using 
0.5588 for the linear term as obtained from Claret et al. 
(2000). The resulting best-fit parameters along with 


iihttps: //www.univie.ac.at/tops/Period04/ 
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Fig. 11 .— Power spectral density plots for the time-series data of four transit residuals. Frequencies have been plotted 
up to the maximum Nyquist frequency. No obvious 1/frequency dependency is visible except (with some good will) for the 
transit of HATP25. 


uncertainties are listed in Table 5. 

For HATP22 we judged that the amount of corre¬ 
lated red noise is minimal. We therefore used TASK8 
(Monte-Carlo) to find parameter uncertainties based on 
5000 simulations. The photometric precision is of high 
quality allowing us to fit for the two limb-darkening 
coefficients simultaneously without producing unphys¬ 
ical results. We also carried out fits where one of the 
parameters were held fixed. All three cases were con¬ 
sistent within the quoted uncertainties. The timing un¬ 
certainty for this light-curve is the smallest with around 
±17 seconds. A close inspection of the HATP22 light- 
curve shows some systematic brightening at around 
BJD 2,456,671.215 and could be due to a star spot 
(Nutzman et al. 2011; Oshagh et al. 2013). However, 
this feature is also present in all the adopted compar¬ 
ison stars and hence is attributable to a short-term 
atmospheric change likely due to thin clouds. 

For HATP12 we again inferred parameter uncertain¬ 
ties using TASK9. Also in this case we judged the er¬ 
rors to be dominated by correlated noise - especially 
prior to planetary ingress. We used a quadratic limb- 


darkening law with initial guess obtained from Claret 
et al. (2000) and listed in Table 4. Both parameters 
were kept freely adjustable during the initial fitting pro¬ 
cess and we did not encounter unphysical results. The 
final best-fit parameters along with their uncertainties 
are shown in Table 5. 

8. Conclusion 

We have presented first results from follow-up ob¬ 
servations of several transiting planets using the Cb- 
NUOJ 0.6m telescope. We have recorded light-curves 
using infocus as well as defocus telescope settings. We 
have qualitatively described various photometric noise 
sources and their effects when operating the telescope 
in a defocus mode. In particular we have outlined the 
advantage of defocus observations and demonstrated to 
achieve a precision of RMS ~ 1 milli-magnitude over 
a time-scale of several hours. We believe the precision 
can be explained by three factors: 1) long integration 
time (enabling the collection of many photons), 2) con¬ 
trolling systematic noise source due to good pointing 
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Table 5. 

Final best-fit parameters for light-curves of HATP25, HATP09, HATP22 and HATP12. The mid-transit time (To) is offset 
by BJD 2,456,000. The uncertainties have been obtained from TASK8 (Monte-Carlo) and TASK9 (Residual-Permutation) 

while setting the adjustable parameter integers to initially 1. 


Parameter 

HATP25 

HATP25 

HATP09 

HATP22 

HATP12 

f a +r b 

o 136 + 0.009 

u.iou_o oio 

0 14 + 0 . 0 i 0 
u * lfl - 0.020 

0 126+ u - uuy 
U.lZU-o.oio 

q 34+U.02 
u -°^-0.02 

0-28++ 

k 

n 197 + 0.003 

u. 10 / -0.003 

PI 1 Q7+0.005 

U.10/-0.004 

0-106™ 

0 106 +0 002 
D. 1 UU- 0.002 

0-1411+0“ 

i(°) 

85 . 3 + q 7 

85.0+ 8 

87+9 

84+ 

88 + 

T 0 (BJD) 

247.12881°;°°°? 

247.129l°;°°J 7 

28 ^ 1671 + 0 - 0005 

ZOO. ID / l_o.oo05 

671.25041°;°°°! 

726 171 n + 0 0004 

/ ZD. 1 / -LU— 0.0004 

U a 

0.3879 (fixed) 

0.3879 (fixed) 

0.5588 (fixed) 

0.61+ 

0 7+0.3 
u - 1 -0.3 

V a 

0.2906 (fixed) 

0.2906 (fixed) 

— 

-o.H8:l 

0 o+ 0 - 7 
u ’ u - 0.8 

ta 

0 1 1 Q+0.008 

u.±±y_ 0 .oo9 

0 191 +0 009 
U. 1 Z 1 -Q 014 

0 114+ 0 - 008 

U.11^1—0.009 

Q O-j + 0.01 

u - 01 -o.oi 

n 24+ 0 - 02 

u - z ^-0.009 

n 

0 016 +0 001 
U.UlU-o.ooi 

0 016 +0 002 
U.UIO —0 002 

0 012 +0 ' 001 
u.uiz_o 001 

0 033+ 0 - 002 
U.UOO —0 002 

n 034+ 0 004 
U.Uolt —0 0007 

a (mmag) 

2.30 

2.30 

1.68 

0.70 

1.85 

Xr 

1.01 

1.01 

1.42 

1.21 

0.91 

JKTEBOP 






uncertainty method 

TASK 8 

TASK9 

TASK 8 

TASK 8 

TASK 8 


precision by stable telescope tracking (keeping the tar¬ 
get star on roughly the same pixels) and 3) beating of 
scintillation noise (as a side-effect from long exposure 
times) which otherwise would be significant for a 0.6m 
telescope using exposure times of a few tenth of seconds 
when operating in in-focus mode. We demonstrated 
that the telescope is well suited for follow-up observa¬ 
tions of suspected transit signals from discovery sur¬ 
veys such as Qatar, HAT, WASP and similar on-going 
projects. However, target stars with possible planets 
have to be bright with a limiting magnitude of around 
V = 13 - 14 (depending also on the transit depth). For 
bright host stars with V ~ 10 a remarkable photometric 
precision of RMS ~ 0.70 milli-magnitude was demon¬ 
strated for HATP22. This precision is high enough 
to determine physical properties of the planet with er¬ 
rors of a few percent. Unfortunately we recorded only 
one complete transit of HATP22. Additional datasets 
would be necessary for a complete analysis also help¬ 
ing to reduce modelling errors. Furthermore, track¬ 
ing of star spots visible during a transit event can be 
used to measure stellar rotation of the host star and 
would be detectable with a sub-milli-magnitude preci¬ 
sion. Mid-transit times were measured to ~ 17 seconds 
for HATP22 to ~ 86 seconds for HATP25 enabling 
some or limited opportunity for transit-timing varia¬ 
tion follow-up studies for the detection of additional 
bodies in the system. In the future we plan to carry 
out follow-up observations of bright host-stars known 
to harbour a transiting planet with the aim to obtain 
~ 1 milli-magnitude photometric precision. 

ACKNOWLEDGMENTS 

The data acquisition and analysis was partially sup¬ 
ported by Basic Science Research Program through the 
National Research Foundation of Korea (NRF) funded 


by the Ministry of Education, Science and Technology 

(2011-0014954). C-H K was partially supported by the 

Korea Research Foundation (KRF) grant funded by the 

Korea Government (NRF-2012R1A1A4A01012467). 

REFERENCES 

Bakos GA, Hartman J, Torres G, Latham DW, Kovacs 
G et ah, ”HAT-P-20b - HAT-P-23b: Four mas¬ 
sive transiting extrasolar planets”, ApJ 742, 116-135 
( 2011 ) 

Burke CJ, McCullough PR, Valenti JA, Long D, Johns- 
Krull CM et al., ”XO-5b: A transiting Jupiter-sized 
planet with a 4 day period”, ApJ 686, 1331-1340 
(2008) 

Charbonneau D, Brown TM, Latham DW, Mayor M, 
’’Detection of Planetary Transits Across a Sun-like 
Star”, ApJ 529, L45-L48 (2000) 

Claret A, ”A new non-linear limb-darkening law for 
LTE stellar atmosphere models. Calculations for - 
5.0 < log[M/H] < +1, 2000 K < T eff < 50000 K 
at several surface gravities”, A&A 363, 1081-1190 
( 2000 ) 

Eastman J, Sivercl R, Gaudi BS, ’’Achieving Bet¬ 
ter Than 1 Minute Accuracy in the Heliocentric 
and Barycentric Julian Dates”, PASP 122, 935-946 
( 2010 ) 

Everett ME & Howell SB, ”A technique for ultrahigh- 
precision CCD photometry”, PASP 113, 1428-1435 
( 2001 ) 

Gillon M, Smalley B, Hebb L, Anderson DR, Triaud 
AHMJ et al., ’’Improved parameters for the tran¬ 
siting hot Jupiters WASP-4b and WASP-5b”, A&A 
496, 259 - 267 (2009) 



12 


T. C. Hinse et al. 


Hartman JD, Stanek KZ, Gaudi BS, Holman MJ, 
McLeod BA, ’’Pushing the Limits of Ground-based 
Photometric Precision: Submillimagnitude Time- 
Series Photometry of the Open Cluster NGC 6791”, 
AJ 130, 2241-2251 (2005) 

Hartman JD, Bakos GA, Torres G, Kovacs G, Noyes 
RW et ah, ”HAT-P-12b: A low-density sub-saturn 
mass planet transiting a nretal-poor K dwarf’, ApJ 
706, 785-796 (2009) 

Henry GW, Marcy GW, Butler RP, Vogt SS, ”A Tran¬ 
siting 51 Peg-like Planet”, ApJ 529, L41-L44 (2000) 

Johns-Krull CM, McCullough PR, Burke CJ, Valenti 
JA, Janes KA et ah, ”XO-3b: A massive planet in 
an eccentric orbit transiting an F5V star”, ApJ 677, 
657-670 (2008) 

Kim CH, Song MH, Yoon JN, Han W, Jeong MJ, ”BD 
Andromedae: A New Short-period RS CVn Eclips¬ 
ing Binary Star with a Distant Tertiary Body in a 
Highly Eccentric Orbit”, ApJ 788, 134-156 (2014) 

Lee JW, Youn J-H, Kim S-L, Lee C-U, Hinse TC, 
’’The Sub-Saturn Mass Transiting Planet HAT-P- 
12b”, AJ 143, 95-100 (2012) 

Lenz P & Breger M, ”Period04 User Guide”, CoAst 
146, 53-136 (2005) 

Mayor M & Queloz D, ”A Jupiter-mass companion to 
a solar-type star”, Nature 378, 355-359 (1995) 

McCullough PR, Burke CJ, Valenti JA, Long D, Johns- 
Krull CM et ah, ”XO-4b: An extrasolar planet tran¬ 
siting an F5V star”, arXiv:0805.2921 (2008) 

Nutzman PA, Fabrycky DC & Fortney JJ, ’’Using Star 
Spots to Measure the Spin-orbit Alignment of Tran¬ 
siting Planets”, ApJ 740, L10-L15 (2011) 

Oshagh M, Santos NC, Boisse I, Boue G, Montalto, 
M et al., ’’Effect of stellar spots on high-precision 
transit light-curve”, A&A 556, 19-24 (2013) 

Quinn SN, Bakos GA, Hartman J, Torres G, Kovacs 
G et ah, ”HAT-P-25b: A hot-Jupiter transiting a 
moderately faint G star”, ApJ 745, 80-89 (2012) 

Shporer A, Bakos GA, Bouchy F, Pont F, Kovacs G 
et al., ”HAT-P-9b: A Low Density Planet Transit¬ 
ing a Moderately Faint F star”, ApJ 690, 1393-1400 
(2009) 

Southworth J, ’’Homogeneous studies of transiting ex¬ 
trasolar planets - I. Light-curve analyses”, MNRAS 
386, 1644-1666 (2008) 

Southworth J, Hinse TC, Jprgensen UG, Dominik M, 
Ricci D et ah, ’’High-precision photometry by tele¬ 
scope defocusing -1. The transiting planetary system 
WASP-5”, MNRAS 396, 1023-1031 (2009) ' 

Young AT, ’’Photometric error analysis. VI. Confirma¬ 
tion of Reiger’s theory of scintillation”, AJ 72, 747- 
744 (1967) 



